Proton spin in a light-cone chiral quark model 
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We discuss the spin structure of the proton in a light-cone treatment of the chiral quark model. 
Based on the fact that the quark helicity (Aq) measured in polarized deep inelastic scattering 
experiments is actually the quark spin defined in the light-cone formalism rather than the quark 
spin (A(/qm) defined in the conventional quark model (or in the rest frame of the nucleon), we 
calculate the x dependence of the polarized quark distribution functions Aq(x), and the polarized 
structure functions gi(x). Special attention is focused on the Melosh-Wigner rotation due to the 
transversal motions of quarks inside the nucleon and Melosh-Wigner rotation effects on the bare 
quark input. It is shown that our results match the experimental data well. 



I. INTRODUCTION 
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It is well known that the nucleon consists of quarks and gluons, and the underlying theory describing their interaction 
is quantum chromodynamics (QCD). However, due to the nonlinear and nonperturbative nature of QCD, it is still 
beyond our practical abilities to calculate from the QCD Lagrangian directly to study the low-energy properties of 
the nucleon at present, and consequently there are many problems that remain unsolved. Among these problems, 
the spin structure of the nucleon is a particularly important one and has triggered a vast number of theoretical and 
experimental investigations. 

In the SU(6) quark model, the proton wave function is given by 
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One can count the number of quarks with certain helicity, =5/3, = — 1/3, eft = 1/3, Gr = 
Naively, one may expect that the sum of quark helicity to be 
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which is equal to the proton spin. However, it was observed in the polarized deep inelastic scattering (pDIS) experiment 
by the European Muon Collaboration |l[ that the sum of quark helicity was much smaller than 1 /2. This is the so-called 
"proton spin crisis" or "proton spin puzzle." Although many models have been proposed to explain this phenomenon, 
the issue is still far from being settled. In particular, we still need theoretical calculation of the x-dependent polarized 
quark distribution functions, with the expectation that these will match the experimental data well, since in most 
previous papers only integrated quantities were discussed, or only parametrizations were given by fitting the data. 

As emphasized by one of us @, the quark spin is not Lorentz invariant, and the quark helicity (Aq) observed 
in pDIS should be viewed as the quark spin defined in the light-cone formalism, which is different from the quark 
spin (Aqqm) defined in the conventional quark model (or in the rest frame of the nucleon). These two quantities are 
connected by the relativistic effect of the Melosh-Wigner rotation 0, [f| due to the transversal motions of quarks inside 
the nucleon. The quantity AS/2 = (Au + Ad + As)/2 represents, in a strict sense, the sum of the quark helicity 
in the light-cone frame rather than the vector sum of the spin carried by quarks and antiquarks in the proton rest 
frame. Hence, the best approach to analyze the experimental data is based on the framework of the light-cone field 



* Corresponding author. Email address: mabqOpku.edu.cn 



2 



theory [6] , or equivalently the infinite momentum technique proposed by Weinberg [7| ■ If we adopt the quantum field 
theory in the standard instant form, the spin of the proton should be the sum of the Melosh-Wigner rotated spin of 
the quarks rather than simply the sum of the spin of the individual constituent quarks. We would like to emphasize 
that although we shall pick a particular model to calculate the quark distribution functions, the relation between Ag 
and AqQM is general, since this effect is a fundamental kinematic relation that must be taken into account regardless 
of the dynamical details [1, HJ. However, not enough attention to this effect was paid in some previous analyses of the 
proton spin. 

In this paper, we derive the polarized quark distribution functions Aq(x) and the polarized structure functions 
<7i (x) for the proton, the neutron, and the deuteron, respectively, in the framework of the chiral quark model without 
referring to polarized experimental data beforehand, and we show that the results can match the experimental data 
well through the fully rclativistic light-cone treatment. The paper is organized as follows. In Sec. [TTl we analyze the 
proton spin in the framework of the chiral quark model. In Sec. IIIH we discuss the bare quark input in the chiral 
quark model and explain the important effect of the Melosh-Wigner rotation. Then we obtain the polarized quark 
distribution functions and the polarized structure functions for the proton, the neutron and the deuteron. In Sec. IIV| 
we summarize the paper. 



II. CHIRAL QUARK MODEL 
A. The basic structure of the chiral quark model 

The chiral quark model [l(l[ll| is widely recognized as an effective theory of QCD at the low-energy scale. It adopts 
a description of its important degrees of freedom in terms of quarks, gluons, and Goldstone bosons at momentum 
scales relating to hadron structure. In this model, the minor effects of the internal gluons are neglected, and the 
valence quarks contained in the nucleon fluctuate into quarks plus Goldstone bosons, which spontaneously break the 
chiral symmetry. This model is successful in explaining a number of problems, including the violation of the Gottfried 
sum rule from the aspect of the flavor asymmetry in the nucleon sea [L2H14| , the NuTeV anomaly resulting from the 
strange-antistrange asymmetry [l5)|. and the isospin symmetry breaking between the proton and the neutron (l6| . It 
is also proposed by Cheng and Li [17j | that the proton spin problem can be accounted for by the quark splitting into 
a quark plus a Goldstone boson in the chiral quark model from an intuitive argument. 

The effective Lagrangian of the interaction between Goldstone bosons and quarks in the leading order is 

C mt = -^7 M 75C> M IT/>, (3) 



where II takes the form 
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Based on the perturbative technique in the light-cone field theory, all particles are on mass shell and the factorization of 
the subprocess is applicable, so we can express the quark distributions inside a nucleon as a convolution of a constituent- 
quark distribution in a nucleon and the structure functions of a constituent quark. The Fock decompositions of 
constituent-quark wave functions are 

\U) = Vz\u )+a w \d-K+) + ^\u7T°)+a K \sK+} + ^\ur } ), 

V2 v 6 

\D) = Vz\d a ) + 0*1*0 + ^\d7T°) + a K \sK°) + %\*q), (5) 

where Z = 1 — §(-Ptt) — (Pk) ~ ^(-P>;) i s t ne renormalization constant for the bare constituent, and |a a | 2 (a = tt, K, rj) 
are the probabilities of finding Goldstone bosons in the dressed constituent-quark states \U) and \D). The fluctuation 
of a bare constituent quark into a Goldstone boson and a recoil bare constituent quark is given as [lH 

j x y v y / 

= Pi^ja ® A& (#)- ( 6 ) 
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Pi^ja(y) is the splitting function, which gives the probability of finding a constituent quark j carrying the light-cone 
momentum fraction y together with a spectator Goldstone boson a, 

where m,, mj, and m a are the masses of the i- and j-constituent quarks and the pseudoscalar meson a respectively. 

Mj a — " J ^" fc ^ + m ^y ± is the square of the light-cone invariant mass of the final state. 

In the effective field theory, it is necessary to introduce an ultraviolet cutoff to make the results finite. The 
conventional choice is to specify the momentum cutoff of the vertex as 
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In this paper, we adopt g' A = 1, following the large N c argument [T^ |. The cutoff parameter A can be determined by 
the experimental data of the Gottfried sum and the constituent-quark-mass inputs for the pion, 
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= \{Z-\{P,) + {Pk) + \{P v )) 
= ^(1 -2 (P.)) 

= 0.235 ± 0.026 (experimental data), (9) 

where (P a ) = Pi^j a (y)dy. We choose the parameters of masses as m u = m d = 330 MeV, m s = 480 MeV, 
= 140 MeV, m K = 495 MeV, m v = 495 MeV, and therefore the cutoff A = 1500 MeV. 



B. Polarized quark distributions in the chiral quark model 

The polarized quark distributions are obtained by replacing the unpolarized splitting function Pi^j a (y) with the 
polarized splitting function APi^j a (y) = P^^j^ a — Pii^jia- The polarized fluctuation process is expressed as 

Aq 3 {x) = f ^-APi^ ja {y)Aq i0 (-) 
Jx V \VJ 

= AP^ a 0A?jo(s). (10) 
If APi^j a (y) is known, one can write down the polarized parton distribution functions 



Au(x) = ZAu {x) + -AP u ^ U7I o ® Au (x) + AP^- <g) Ad (x) + AP U ^ UV <g> Au (x), 



Ad(x) = ZAd a (x) + -AP d ^ d7r o ® Ad (x) + AP u _> d7T + <g> Au (x) + AP d ^ dv ® Ad (x), 
As(x) = AP u ^ sK + ® Au„{x) + AP d ^ sK o ® Ad (x), 

Au(x) = Ad{x) = As(x) = 0. (11) 

We should point out that in the chiral quark model, antiquarks are only produced in the process of the splitting 
of Goldstone bosons. Since Goldstone bosons are unpolarized, the polarization of antiquarks is zero, and this is 
compatible with the available experimental data (20l |. 

In most models including the chiral quark model, the possibility of "helicity reverse" is correlated with the matrix 
element 

~7^ 15 ^/^' (12) 
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FIG. 1: The behavior of the unpolarized and polarized splitting functions P and AP versus x. 



which describes the change of the helicity of a fermion after interacting with a pseudoscalar particle. It is often taken 
as an assumption [17] that the helicity of the fermion is totally reversed, and therefore APi^>j a (y) = —Pi^>j a (y) in 
the explanation of the proton spin problem. However, as stressed in the introduction, we should adopt the light-front 
formula @ to analyze the spin structure of a composed system carefully. In this approach we can obtain 

u{k,m,\) u{p,M,\) 2 r , , /,! . 2W i . . 2M no\ 

^+ 75 -j= = ^+ [-k+p+ + mM - (k -ik){p +ip)\, (13) 

u(k,m,t) u (p, M, j) 2 r . x . 2 , , ! 

75 = -^T^T L m (P - *P ) + M ( fc - * fc )J > ( 14 ) 

^+ 75 ^= = - mM +(k +ik ){p -ip)\, (15) 

u(k,m,i) w(p,M,f) 2 r , ! . 2 s - , 2M rift s 

75 = -^+^+ L m (P + *P ) + M ( k +** )J • (16) 

Furthermore, the following kinematics is generally associated with the vertex 

P=(^,0x), fc=L + ,!l±^ ffc i|. (17) 
We obtain the form of the helicity reverse function APr, 

(—VP +2 + m 4 m,) 2 — m 2 k~± 2 , v 

AP* = P (t^t - t=H0 = f-^- ^ i4r, (18) 

{—yp + + rriirrij) + mfk± 2 

which is not — 1 . Accordingly, in the chiral quark model, the polarized splitting function APj_>.j a (y) takes the form 

AP^_ yj a (y) = Pi— >j Q 

1 f gAjm +m j ) \ 2 f a (rrij - m^) 2 - kj_ 
- 32ttH / ; 7 X J/ 2 (l-y)[mf-M? Q p- liJj 

The behaviors of P and AP are displayed in Fig. [TJ 

It is important to distinguish two different types of quark and gluon contributions to the nucleon sea measured in 
experiments. One is the perturbative "extrinsic" sea, which is generated from the QCD hard bremsstrahlung and 
the gluon-splitting process, and this type of sea is associated with the internal composition of gluons rather than 
the proton itself. The other one is the nonperturbative "intrinsic" sea, which is directly connected to the valence 
quarks of the nucleon. The quark distributions measured at certain Q 2 include not only the intrinsic sea but also 
the extrinsic sea (2lj . However, in the chiral quark model, only the intrinsic sea is calculated while the extrinsic sea 
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is left out. Therefore, it is generally necessary to add the extrinsic sea before comparing results with experiments 
when we consider the unpolarized quark distributions [l4j ]. Nevertheless, when we consider the polarized quark 
distributions, we assume that the perturbative extrinsic sea is totally unpolarized, since the correlative processes give 
equal contributions to both types of polarization. Therefore we do not need to add the perturbative extrinsic sea. 



III. THE BARE QUARK INPUT AND THE MELOSH-WIGNER ROTATION 

A. The bare quark input 

We need the input of polarized bare quark distributions Auq(x) and Ado(x) in the calculation. However, since 
these distributions are not measurable in experiments, there is no parametrization available. It is also not proper 
to use polarized parton distribution parametrization of valance quark as the input directly, since these two kinds of 
distributions are physically different. In this paper, we try to determine them from the unpolarized quark distributions, 
since the chiral quark model should give a unified description of the unpolarized and the polarized quark distributions. 
Therefore, we actually use the chiral quark model as a theoretical tool to connect the unpolarized and the polarized 
quark distributions, and we never refer to polarized experimental data beforehand. 

We assume that the unpolarized bare quark distributions uq(x) and do(x) take the form 

_ 2x a "(l~x) K 
Mx) BK + 1A + I]' 
, , n x a "(l~x) bd 

B [ad + l,b d + l\ 

where B[i,j] is the Euler beta function. It is observed that q(x) — > Zqo(x) when x — > 1 in the chiral quark model, so 
we can determine b u and bd from the asymptotic relation given by the unpolarized parametrization. Afterward, we 
determine the value of a u and ad by requiring that the maximum points of functions xq{x) {q = u, d) should be at 
the same values of x as those in the parametrization. It is found that the approximate values for these parameters 
are a u = 0.5, ad = 1, b u = 3, b u = 5 based on the Martin-Roberts-Stirling-Thorne parametrization [22J at fixed 
Q 2 = 2.5 GeV 2 . Then we can adopt the naive SU(6) quark model to get the spin structure of bare quarks in the 
proton rest frame, 

Au$ M (x) = ^u (x), Ad$ M (x) = -~d (x). (21) 

Before obtaining the result in the light-cone frame, we need to analyze the effect of the Melosh-Wigner rotation 
carefully. 



B. The Melosh-Wigner rotation 

The Melosh-Wigner rotation [J, |5J is a natural result of relativistic kinematics, and was shown to play an essential 
role in the analysis of the spin structure @, Q . The key point lies in the fact that the vector sum of the constituent 
spin for a composite system is not Lorentz invariant from the relativistic viewpoint. The quark helicity distribution 
function Aq measured in pDIS is defined by the matrix element 

Ag=(p t |?7 + 75gb t }- (22) 

If we express the quark wave functions in terms of light-cone Dirac spinors @ , we will obtain 

Aq(x) = qT(x)-q±(x), (23) 

where q^(x) and q^(x) are probabilities of finding a quark of flavor q with fraction x of the proton longitudinal 
momentum and with polarization parallel and antiparallel to the direction of the proton spin in the proton infinite 
momentum frame (i.e., the quark spin states in the light-cone frame). However, if we express the quark wave functions 
in terms of conventional instant form Dirac spinors (i.e., the quark spin states in the proton rest frame), we will get 
that 



Aq(x) = 




J d 2 k~lw D (x, kT) A9QA/ (x, kT) , (24) 
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with 



/ — - \ (k + + m) 2 — I: - 
W D (x,k ± )=i L 4 (25) 

V 7 (fc+ + m) 2 + k ± 2 

being the contribution from the rclativistic effect due to quark transversal motions, q Sz= L (x, k±j and ? Ss , = _ i (x, k±\ 

being probabilities of finding, in the proton rest frame, a quark or antiquark of flavor q with rest mass m and 
momentum p^ and with spin parallel and antiparallel to the rest proton spin respectively, and k + = xM. with 

M 2 = J2 t (m 2 + hl 2 )/x i . 

It is straightforward to express the quark spin in the light-cone frame with the effect of the Melosh-Wigner rotation 

as 

Au {x) = Au% M (x)W u {x) = ~u (x)W u (x), 

Ad Q {x) = Ad$ M {x)W d {x) = ~d {x)W d {x), (26) 

where W u (x) and W d (x) are Melosh-Wigner rotation factors for the u quark and the d quark, respectively. In this 
paper, we adopt m q = 330 MeV, m Du = 600 MeV, m Dd = 900 MeV, a u = 330 MeV, and a d = 200 MeV as inputs 
for the Melosh-Wigner rotation factor. The Melosh-Wigner rotation factors for the u quark and the d quark are 
displayed in Fig. [3] The difference of inputs for the Melosh-Wigner rotation between u quark and d quark results 
from the difference between the u quark and the d quark in the wave function ([1]). When the u quark is probed, the 
spectator quarks have both the scalar (S — 0) and the vector (5* = 1) components, but when the d quark is probed, 
the spectator quarks only have the vector (5* = 1) component. Therefore, the u quark has a more significant effect 
than the d quark under Melosh-Wigner rotation. 



C. The polarized quark distribution functions 



With all the above considerations taken into account, we can get the polarized quark distribution functions Au(x), 
Ad(x), and As(x), and the results are displayed in Figs. [31 HI El It should be pointed out that the experimental data 
have been evaluated to a common Q 2 — 2.5 GeV 2 in order to compare with the theoretical prediction, which is also 
conducted at the same Q 2 . It can be found that our results match with the experimental data well. 

Furthermore, we display the polarized structure functions g± ,(N — p, n, d) in Figs. El M It is shown that when 
we use the chiral quark model with the effect of the Melosh-Wigner rotation taken into consideration, the behaviors of 
these functions can match the experimental data well. In contrast, when the Melosh-Wigner rotation is not included, 
the description of the behavior in large x region is not satisfying. Hence, the Melosh-Wigner rotation can provide an 
important mechanism to depress the polarization in large x region. 

To make our analysis more comprehensive and persuasive, we also evaluate the polarized quark distribution func- 
tions to a much higher energy scale Q 2 = 50 GeV 2 with the Higher Order Perturbative Parton Evolution Toolkit 
(HOPPET) [13, and show the behavior of at Q 2 = 50 GeV 2 in Figs. H ® respectively. We find that the 
Q 2 -evoluted results can match better with the experimental data at large x, in which region Q 2 is much larger than 
2.5 GeV 2 , and also the results are insensitive to different inputs of gluon contributions. However, the experimental pre- 
cision is still not good enough to distinguish the effect of Q 2 evolution and we expect more experimental investigations 
concerning the QCD evolution behaviors of spin related quantities of the nucleon in the future. 



IV. SUMMARY 



In this paper, we discuss the spin structure of the proton in the relativistic light-cone formulae. We stress that the 
quark helicity measured in pDIS is actually the quark spin defined in the light-cone formalism and is different from 
the quark spin defined in the quark model (or rest frame of the nucleon), and these two quantities are related by 
the relativistic kinematical effect of the Melosh-Wigner rotation. This effect must be taken into account regardless 
of the dynamical details. We adopt the chiral quark model to calculate the polarized quark distribution functions 
Aq(x) and the polarized structure functions g%. We show that the results can match the experimental data well. 
However, we would like to emphasize that the chiral quark model may possibly be replaced by some other kinds of 
dynamics, but the relativistic kinematics of the Melosh-Wigner rotation is prerequisite. We suggest that more precise 
experiments should be carried out to enable more accurate determination of polarized quark distributions to confront 
with theoretical calculations. 
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FIG. 4: The polarized quark distribution function for d quark xAd(x) versus x. The solid line is the result from the chiral quark 
model with Melosh-Wigner rotation taken into account. The dashed line is the result from the pure chiral quark model without 
taking Melosh-Wigner rotation into account. The dotted line is the result from the chiral quark model using the nonrelativistic 
vertex. The experimental data of HERMES are from Ref. [2(| and evaluated at a common Q 2 = 2.5 GeV 2 . 
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FIG. 5: The polarized quark distribution function for s quark xAs(x) versus x. The solid line is the result from the chiral quark 
model with Melosh-Wigner rotation taken into account. The dashed line is the result from the pure chiral quark model without 
taking Melosh-Wigner rotation into account. The dotted line is the result from the chiral quark model using the nonrelativistic 
vertex. The experimental data of HERMES are from Ref. [2(| and evaluated at a common Q 2 = 2.5 GeV 2 . 




FIG. 6: The polarized structure function for the proton g\ versus x. The solid line is the result from the chiral quark model 
with Melosh-Wigner rotation taken into account. The dashed line is the result from the pure chiral quark model without taking 
Melosh-Wigner rotation into account. The dotted line is the result from the chiral quark model using the nonrelativistic vertex. 
The dash-dotted line is the result after evolution to Q 2 = 50 GeV 2 . The experimental data of HERMES, E143, and COMPASS 
are from Refs. [2ul.E^l 25 . respectively. The left and right figures differ in the x-axis type for the purpose of showing the small- a; 
(left figure) and large- a; (right figure) behaviors. We only include data with x > 0.1 in the right figure. 
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FIG. 7: The polarized structure function for the neutron g" versus x. The solid line is the result from the chiral quark model 
with Melosh-Wigner rotation taken into account. The dashed line is the result from the pure chiral quark model without taking 
Melosh-Wigner rotation into account. The dotted line is the result from the chiral quark model using the nonrelativistic vertex. 
The dash-dotted line is the result after evolution to Q 2 = 50 GeV 2 . The experimental data of HERMES and E143 are from 
Refs. H3,[2^], respectively. The left and right figures differ in the a;-axis type for the purpose of showing the small- a; (left figure) 
and large- x (right figure) behaviors. We only include data with x > 0.1 in the right figure. 




FIG. 8: The polarized structure function for the deuteron gf versus x. The solid line is the result from the chiral quark model 
with Melosh-Wigner rotation taken into account. The dashed line is the result from the pure chiral quark model without 
taking Melosh-Wigner rotation into account. The dotted line is the result from the chiral quark model using the nonrelativistic 
vertex. The dash-dotted line is the result after evolution to Q 2 = 50 GeV 2 . The experimental data of HERMES, E143, SMC, 
and COMPASS are from Refs. [2p|. l23l - [26l ] . respectively. The left and right figures differ in the x-axis type for the purpose of 
showing the small- a; (left figure) and large- a; (right figure) behaviors. We only include data with x > 0.1 in the right figure. 



